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1,0  INTRODUCTION 

This  Annual  Scientific  Report  sumnarizes  research  accomplished  during 
the  first  year  (January  82  to  January  83)  of  AFOSR  Grant  No.  AF0SR-82-0041 , 
The  objective  of  this  effort  is  to  determine  the  manner  in  which  "small" 
fatigue  cracks  initiate  at  notches,  extend  by  cyclic  loading,  Interact 
with  adjacent  flaws,  and  coalesce  Into  a  single  dominant  crack  which  con¬ 
trols  final  fracture,  The  desired  product  is  a  predictive  scheme  capable 
of  analyzing  the  early  stages  of  fatigue  crack  growth  which  are  character¬ 
ized  by  the  growth  and  link-up  of  small  cracks.  Research  toward  this 
goal  is  directed  at  the  following  tasks. 

I.  Crack  Growth  Predictive  Algorithm.  A  computer  program  will  be  deve¬ 
loped  to  predict  the  growth  and  coalescence  of  multiple  cracks  loca¬ 
ted  at  notches, 

II.  Crack  Interaction  Analysis.  Stress  intensity  factors  solutions  will 
be  computed  for  multiple  cracks  located  at  an  open  hole.  These 
solutions  are  required  for  the  multiple  crack  growth  algorithm. 

III.  Crack  Coalescence  Experiments.  Fatigue  tests  will  be  conducted  with 
multiply  cracked  specimens  to  provide  a  data  base  to  evaluate  the 
predictive  model.  Initially  the  model  will  be  verified  with  "large 
crack"  results  directed  toward  coalescence  aspects  of  the  problem. 
Subsequent  experiments  will  focus  on  coalescence  of  "small"  cracks. 

IV.  Characterization  of  Small  Cracks.  This  phase  of  the  effort  Is 
directed  toward  the  growth  and  coalescence  of  physically  small  cracks. 
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These  date  will  be  used  to  characterize  the  growth  of  "small"  flaws, 
which  are  expected  to  behave  differently  than  "large"  cracks. 

Additional  details  of  these  goals  and  progress  to  date  are  described  in 
the  remaining  sections  of  this  report. 


2,0  SUMMARY  OF  PROGRESS  AND  ACCOMPLISHMENTS 

2.1  Task  I  Progress  "  Crack  Growth  Algorithm 

A  computer  program  has  been  written  to  analyze  the  multiple  crack 
configurations  shown  in  Fig,  1,  Various  combinations  of  surface  and/or 
corner  flaws  are  assumed  to  be  located  along  the  bore  of  a  circular  hole. 
It  is  assumed  that  linear  elastic  fracture  mechanics  are  valid  and  that 
crack  growth  rates  are  controlled  by  the  stress  Intensity  factor  K. 

Stress  intensity  factors  are  estimated  at  the  major  and  minor  crack  axes 
(locations  1  to  6  in  Fig,  1)  by  modifying  the  solutions  for  single  cracks 
in  a  large  body  (1)  with  an  interaction  factor  (2)  developed  here  for 
coalescing  cracks,  (The  interaction  solution  is  discussed  in  the  Task 
II  progress  section,)  The  crack  tips  are  allowed  to  grow  independently 
at  their  major  and  minor  axes  without  specifying  crack  shape  changes. 

The  multi -degree  of  freedom  program  predicts  the  growth  of  the  separate 
initial  cracks  and  their  coalescence  into  single  surface  corner,  or 
through-the- thickness  flaws  as  shown  in  Fig.  2.  Sample  predictions  and 
comparisons  with  test  data  are  described  in  the  Task  III  progress  section. 

2.2  Task  II  Progress  -  Crack  Interaction  Analysis 

A  key  element  of  the  overall  predictive  scheme  is  determining  the 
interaction  between  adjacent  cracks  prior  to  their  actual  coalescence. 

It  is  expected,  for  example,  that  the  presence  of  a  nearby  flaw  will 
cause  an  increased  crack  growth  rate.  Although  a  few  authors  have 
studied  the  magnification  in  stress  intensity  factors  due  to  coalescing 
Cracks  (3-6),  solutions  were  not  available  for  multiple  cracks  at  notches. 
Thus,  one  phase  of  the  current  effort  is  directed  toward  obtaining  stress 
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intensity  factor  solutions  for  coalescing  cracks  at  fastener  holes,  a 
geometry  which  is  also  studied  experimentally  in  the  current  program. 

The  three-dimensional  finite  element-alternating  method  (7-9)  was 
used  to  compute  stress  Intensity  factors  for  symmetric  corner  cracks 
located  at  opposite  sides  of  a  plate  containing  a  hole  loaded  in  remote 
tension  as  shown  in  Fig,  3,  Although  the  analysis  was  limited  to  sym¬ 
metric  corner  cracks,  several  flaw  shapes  and  sizes  were  considered. 

By  comparing  the  double  crack  results  with  corresponding  solutions  for 
single  cracks,  it  was  possible  to  determine  the  effect  of  the  second 
crack  on  the  stress  intensity  factor  as  the  two  flaws  grow  toward  each 
other  and  eventually  touch,  The  influence  of  the  adjacent  flaw  was 
represented  as  an  interaction  factor  y ,  The  y  factor  is  defined  as  the 
ratio  of  the  double  crack  K  divided  by  the  single  crack  result,  and  is 
given  in  Figs,  4  and  5  for  two  positions  along  the  crack  border. 

Note  in  Fig.  4  that  a  significant  Increase  in  K  [y  >  1)  occurs 
at  the  hole  bore  position  (<j>  =  90°  in  Fig,  3)  as  the  adjacent  cracks 
approach  each  other  (as  the  dimensionless  crack  separation  distance 
t$ep/a  ■*  0)-  This  increased  stress  intensity  factor  is  quite  localized 
along  the  crack  perimeter,  however,  and  as  shown  In  Fig.  5,  is  much  less 
significant  at  the  front  surface  location  ($  =  0  in  Fig,  3).  Note  in 
Fig,  5,  for  example,  that  there  is  actually  a  slight  reduction  in  K 
(y  <  1)  at  the  free  surface  as  the  flaws  approach  each  other.  Although 
the  interaction  Is  a  strong  function  of  crack  spacing  (as  represented  by 
dimensionless  spacing  tsep/a  In  Figs.  4  and  5),  crack  shape  a/c  has 
relatively  little  effect  for  the  cases  examined  here. 
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Additional  details  of  the  crack  interaction  analysis  are  given  in 
Refs.  2  and  10.  Reference  2  is  included  in  this  report  as  Appendix  A, 

The  interaction  factors  obtained  in  this  Task  are  incorporated  in  the 
computer  program  described  in  Task  1  to  account  for  crack  coalescence 
in  the  overall  prediction  scheme.  Experimental  verification  of  the 
crack  coalescence  predictions  are  described  in  the  following  section. 

2.3  Task  III  Progress  -  Crack  Coalescence  Experiments 

This  phase  of  the  effort  Is  aimed  at  evaluating  the  numerical  scheme 
described  in  Task  I  for  predicting  the  initial  growth  and  coalescence  of 
multiple  cracks  located  at  notches.  In  order  to  separate  crack  coalescence 
aspects  from  "small  crack"  phenomena,  the  initial  work  has  Involved 
experiments  with  coalescence  of  relatively  "large"  cracks.  Current  work 
deals  with  the  small  crack  problem  and  is  described  separately  in  the 
Task  IV  section, 

2.3,1  Single  Crack  Verification 

Prior  to  describing  crack  coalescence  experiments,  Initial  verifi¬ 
cation  of  the  predictive  model  for  single  flaws  will  be  discussed.  In  an 
earlier  program,  fatigue  crack  growth  tests  were  conducted  with  large 
specimens  machined  from  polymethylmethacrylate  (PMMA),  a  transparent 
polymer.  The  test  specimen,  mounted  in  loading  grips,  is  shown  schematic 
cally  In  Fig.  6,  The  specimens  were  0,7  in.  thick,  8,0  In.  wide,  contained 
a  0.75  in,  hole,  and  were  loaded  In  remote  tension,  Semi elliptical  sur¬ 
face  cracks  or  quarter  elliptical  corner  cracks  were  located  along  the 
hole  bore  In  the  configuration  shown  in  Fig,  2,  A  camera  viewed  the 
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crack  plane  through  the  transparent  specimen  by  means  of  a  mirror  placed 
in  the  cutout  area  of  the  grips.  Crack  growth  was  photographed  as  a  func¬ 
tion  of  elapsed  cycles,  and  the  photos  were  measured  to  give  the  crack 
dimensions. 

Typical  results  for  a  semielliptical  crack  located  midway  along  t^e 
hole  bore  are  given  in  Fig,  7,  Fig.  7a  presents  crack  growth  prior  „o 
penetration  through  the  specimen  thickness  by  the  hole  bore  crack  c  -nen- 

T 

sion  2a.  (Thickness  penetration  occurs  when  a  *  ^  *  0.35  inch  in  Fig.  7a.) 
Figure  7b  describes  crack  growth  during  the  transition  period  when  the 
part-through  crack  grows  into  a  uniform  through-the-thickness  flaw. 

The  experimental  measurements  are  represented  by  open  symbols  in  Fig.  7, 
while  predictions  based  on  a  single  crack  version  of  the  predictive  pro¬ 
gram  are  given  by  the  solid  lines.  Note  that  the  numerical  algorithm 
predicts  the  crack  growth,  Including  shape  changes,  quite  well.  Additional 
comparisons  with  single  surface  and  corner  crack  tests  are  given  in  Ref,  11. 

Thus,  confident  that  single  cracks  could  be  predicted  fairly  well, 
the  present  program  is  studying  the  coalescence  problem,  The  following 
subtasks  are  being  examined, 

2.3,2  Multiply  Cracked  Holes  In  PMMA  Plates 

Several  fatigue  tests  have  been  conducted  with  PMMA  specimens  modeled 
after  the  double  comer  crack  geometry  studied  by  the  finite  element- 
alternating  method  described  in  Task  II.  The  test  specimens  were  3.5 
Inches  wide,  had  a  0,375  inch  diameter  hole,  and  were  0,7  Inch  thick. 

Two  comer  cracks  were  Introduced  f  opposite  sides  of  the  plate  thickness, 
to  match  the  symmetric  com*  crock  configuration  shown  in  Fig,  3, 
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The  specimens  were  fastened  to  metal  grips  similar  to  those  described 
previously  In  Fig.  6,  and  crack  growth  was  again  photographed  through 
the  transparent  end  of  the  specimen.  Fatigue  crack  growth  curves  are 
shown  for  two  specimens  in  Figs,  8  and  9,  Here  al  and  cl  represent  the 
hole  bore  and  free  surface  dimension  of  one  crack  (see  Fig.  lc),  while 
a2  and  c2  are  the  corresponding  dimensions  for  the  second  flaw.  Note 
that  the  two  cracks  were  quite  symmetric  in  these  two  tests,  as  assumed 
in  the  Task  II  stress  intensity  factor  analysis.  After  crack  coalescence 
(al  +  a2  *  T  =  specimen  thickness),  only  the  free  surface  dimensions 
(cl  and  c2)  are  given  In  Figs.  8  and  9. 

The  solid  and  dashed  lines  give  predictions  by  the  multiple  degree 
of  freedom  model  described  earlier.  The  dashed  lines  predict  the  c 
dimensions  prior  to  crack  coalescence,  while  the  solid  lines  represent 
predictions  for  the  a  dimension.  The  single  solid  line  following  pre¬ 
dicted  coalescence  give  the  estimated  growth  of  a  single  through-the- 
thickness  flaw.  Although  the  predictions  are  somewhat  conservative, 
they  provide  a  reasonable  estimate  for  total  specimen  life.  It  Is 
planned  to  conduct  additional  tests  with  nonsymmetric  corner  cracks, 

2.3.3  Multiple  Cracks  in  PMMA  Bend  Specimens 

Several  crack  coalescence  tests  have  been  conducted  with  the  semi¬ 
circular  notched  bend  geometry  shown  In  Fig,  10,  The  bend  geometry  was 
selected  since  loading  can  be  applied  without  bonding  grips  to  the 
specimen  ends,  and  the  crack  plane  can  be  observed  directly  through 
the  transparent  PMMA  without  use  of  a  viewing  mirror.  Again  corner 
cracks  were  located  on  opposite  sides  of  the  notch  (see  Fig,  10b)  and 


time  lapse  photography  was  used  to  measure  and  record  crack  growth. 

A  typical  set  of  photographs  showing  initial  crack  growth  and 
coalescence  is  given  in  Fig.  11.  Note  that  crack  coalescence  occurs 
at  approximately  11,890  cycles.  The  cusp  dimension  h  (see  Fig.  10c) 
then  grows  very  rapidly  until  a  fairly  uniform  through- the-thickness 
crack  front  is  achieved  at  approximately  12,200  cycles.  Typical  plots 
of  measured  crack  dimensions  are  given  as  functions  of  elapsed  cycles 
in  Figs.  12-15.  It  is  interesting  to  note  that  in  several  experiments 
the  crack  tips,  did  not  coalesce  immediately^  but  grew  past  each  other 
in  slightly  different  planes  to  form  the  overlapped  configuration 
shown  schematically  in  Fig.  16.  Current  efforts  are  directed  toward 
estimating  stress  intensity  factors  for  the  notched  bend  geometry  so 
that  the  life  prediction  model  can  be  applied  to  this  series  of  tests. 

2.3.4  Multiple  Cracks  in  Metal  Edge  Notch  Specimens 

The  multi -degree  of  freedom  fatigue  crack  growth  model  has  also 
been  evaluated  with  experimental  data  provided  by  Dr.  A.  Thakker  (12). 
He  conducted  a  series  of  tests  with  the  double  notched  specimen  shown 
in  Fig.  17.  The  specimens  were  made  from  Waspalloy,  a  turbine  disk 
alloy,  and  were  tested  at  400°F.  Small  initial  cracks  were  introduced 
along  the  bore  of  one  of  the  semicircular  edge  notches,  and  heat 
tinting  techniques  were  used  to  periodically  mark  the  fatigue  crack 
fronts.  Crack  length  and  shape  measurements  were  then  made  on  the 
fracture  surface  following  final  failure. 


The  specimen  dimensions,  load  history,  and  crack  growth  properties 
were  provided  along  with  the  measured  crack  lengths.  The  open  hole  stress 
intensity  factor  solution  employed  in  the  original  crack  growth  model 
was  modified  here  to  reflect  the  fact  that  the  cracks  occured  at  edge 
notches  rather  than  central  holes.  The  multiple  crack  computer  program 
was  then  used  to  predict  the  growth  of  the  individual  cracks,  coalescence, 
and  final  fracture. 

The  results  for  several  specimens  are  given  in  Figs.  18-23.  Note 
that  since  fracture  surface  markings  were  used  to  determine  the  crack 
lengths,  relatively  few  experimental  data  points  are  shown.  In  addition, 
few  details  are  known  regarding  crack  coalescence  and  flaw  shape  changes. 
Nevertheless,  when  the  total  crack  growth  picture  is  viewed,  the  predic¬ 
tive  model  generally  does  an  excellent  job  in  predicting  total  specimen 
life. 


2.4  Task  IV  Progress  -  Characterization  of  Small  Cracks 

It  is  commonly  known  that  physically  small  cracks  (less  than  0.10 
inch  long)  may  grow  faster  than  larger  flaws  subjected  to  similar  stress 
intensity  factor  levels.  The  assumptions  regarding  continuum  behavior  and 
linear  elastic  material  response  may  lose  validity  as  crack  lengths 
become  increasingly  small.  The  problem  is  further  complicated  by  the 
presence  of  adjacent  flaws  which  often  initiate  simultaneously  and 
which  interact  prior  to  the  coalescence  process.  The  preceeding  sections 
have  described  tasks  directed  toward  the  coalescence  problem. 

This  prior  work  has  concentrated  on  relatively  large  flaws  in  order 
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to  separate  small  crack  behavior  from  coalescence  phenomena.  This 
section  briefly  outlines  work  aimed  at  the  "small  crack  problem." 

A  series  of  fatigue  tests  have  been  conducted  with  V-notched  PMMA 
specimens  as  shown  schematically  in  Fig.  24.  Again  4-point  bending 
was  applied  to  simplify  specimen  loading  and  crack  plane  observation. 
Natural  fatigue  cracks  were  allowed  to  develop  naturally  along  the 
V-notch  and  to  coalesce  into  a  single  flaw.  Time  lapse  photography 
was  used  to  record  and  measure  the  individual  and  coalescing  flaws. 

Typical  results  for  two  of  the  crack  initiation  tests  are  given 
in  Figs.  25  and  26.  Note  that  width  and  depth  measurements  are  given 
as  functions  of  elapsed  cycles  for  the  naturally  occurring  cracks. 

Crack  spacing  is  also  measured  (not  shown  in  Figs.  25  and  26)  so  that 
coalescence  behavior  can  be  studied  as  well.  Flaw  sizes  less  than 
0.005  inches  are  routinely  observed  and  measured. 

In  several  tests,  multiple  cracks  formed  uniformly  along  the  notch 
and  quickly  joined  to  form  a  short  through-the- thickness  crack.  Results 
for  two  of  these  experiments  is  given  in  Figs.  27  and  28.  Here  the 
initial  through-the- thickness  cracks  size  is  on  the  order  of  0.003  in. 

Analysis  of  these  small  crack  experiments  is  currently  in  progress. 
Analysis  of  the  notched  specimen  are  being  performed  to  determine  the 
notch  tip  stresses.  It  is  hoped  that  the  through-crack  results  can 
be  used  to  determine  baseline  properties  for  small  crack  behavior  and 
that  the  Task  I  predictive  algorithm  can  be  modified  to  predict  the 
multiple  crack  growth  and  coalescence  results.  It  is  recognized  that 
this  goal  could  be  complicated  by  nonlinear  plastic  zone  effects  which 
could  limit  the  LEFM  approach  for  the  small  crack  tests. 


3.0  CURRENT  AND  PLANNED  RESEARCH 

As  indicated,  current  efforts  are  directed  toward  analysis  of  the 
semicircular  notch  coalescence  experiments  (2.3.3)  and  the  small  crack 
results  (2.4).  It  is  also  planned  to  conduct  additional  tests  with 
cracked  hole  specimens  (2.3.2)  which  contain  nonsymmetric  corner  cracks. 

The  prior  results  have  all  been  obtained  with  PMMA  specimens.  Since 
PMMA  is  a  relatively  brittle  material,  tests  will  also  be  conducted  with 
polycarbonate,  a  more  ductile  polymer,  in  order  to  emphasize  the  plasti¬ 
city  problem.  A  sheet  of  polycarbonate  has  been  purchased  and  several 
preliminary  characterization  tests  have  been  conducted.  It  is  planned 
to  repeat  the  PMMA  coalescence  experiments  with  the  polycarbonate  speci¬ 
mens.  A  modified  V-notch  specimen,  with  a  slightly  larger  notch  radius, 
is  being  considered  to  minimize  variations  in  the  current  "sharp" 

V-notch. 

Following  the  remaining  analysis  of  the  PMMA  results  and  the  poly¬ 
carbonate  tests,  it  was  originally  planned  to  conduct  a  small  study  with 
metal  rpecimens.  Since  the  current  predictive  model  has  matched  existing 
metal  results  quite  well  (2.3.4),  and  since  it  is  much  more  difficult 
to  measure  internal  crack  dimensions  with  metal  specimens,  it  may  be 
desirable  to  continue  the  transparent  polymer  study.  Overload/ retardation 
experiments  could  be  conducted  with  the  polycarbonate  material,  for 
example,  to  determine  load  history  effects  on  the  coalescence  process. 

(All  tests  conducted  to  date  have  employed  constant  amplitude  cyclic 
loads).  Since  the  crack  tip  plastic  zone  varies  through  the  specimen 
thickness,  and  crack  interaction  effects  vary  along  the  flaw  perimeter 
(2.2),  one  would  expect  crack  coalescence  to  be  a  more  complex  three- 
dimensional  problem  in  ductile  materials  subjected  to  variable  load 
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histories.  The  decision  to  conduct  variable  load  tests  will  be  discussed 
with  the  project  monitor  following  the  outcome  of  the  current  PMMA  and 
polycarbonate  tests. 
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4.0  PROFESSIONAL  PERSONNEL,  PUBLICATIONS,  AND  PRESENTATIONS 

4.1  Personnel 

A. F.  Grandt,  Jr.  has  served  as  principal  investigator  for  the  research 
effort.  He  has  beer  assisted  by  the  following  graduate  research  assistants. 

B. J.  Heath  Mr.  Heath  completed  a  Masters  Thesis  entitled  "Stress 

Intensity  Factors  for  Coalescing  and  Single  Corner  Flaws  Along 
A  Hole  Bore  in  a  Plate"  in  January  of  1983  and  was  officially 
awarded  the  Degree  of  Master  of  Science  in  Aeronautics  and 
Astronautics  in  May  of  1983. 

R.  Perez  expects  to  complete  requirements  for  the  Master  of  Science 
Degree  in  Aeronautics  and  Astronautics  in  May  of  1983  with  the 
Thesis  entitled  "Cyclic  Growth  and  Coalescence  of  Multiple 
Fatigue  Cracks  Located  at  Notches." 

J.E.  Pope  Mr.  Pope  is  currently  pursuing  a  Ph.D.  program  in 
Aeronautics  and  Astronautics. 

D.E.  Tritsch,  a  Teaching  Assistant  supported  by  other  School  of 
Aeronautics  and  Astronautics  Department  funds,  has  also  worked 
on  a  Masters  Thesis  which  supports  the  goals  of  this  research 
grant. 

4.2  Publications  and  Presentations 

The  following  reports  have  been  prepared  for  publication 

B.J.  Heath  and  A.F.  Grandt,  Jr.,  "Stress  Intensity  Factors  for 

Coalescing  and  Single  Corner  Flaws  Along  a  Hole  Bore  in  a  Plate," 

Engineering  Fracture  Mechanics, (in  press). 

B.J.  Heath,  "Stress  Intensity  Factors  for  Coalescing  and  Single 

Comer  Flaws  Along  a  Hole  Bore  in  a  Plate,"  M.S.  Thesis,  Purdue 

University,  May  1983. 
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R.  Perez  and  A.F.  Grandt,  Jr.,  "Coalescence  of  Multiple  Fatigue 
Cracks  at  a  Notch,"  Proceedings  of  the  American  Society  of 
Civil  Engineers  Engineering  Mechanics  Division  Speciality  Conference, 
Purdue  University,  May  1983.  Paper  to  be  presented  24  May  1983. 

The  following  publications  are  In  preparation. 

A.F.  Grandt,  Jr.,  R.  Perez,  D.  Trltsch,  and  A.  Thakker,  "Fatigue 
Life  Predictions  for  Multiply  Cracked  Notches,"  planned  for 
submission  to  Journal  of  Engineering  Materials  and  Technology. 

D.E.  Trltsch,  "Fatigue  Crack  Size  and  Shape  Predictions,"  M.S.  Thesis, 
Purdue  University,  expected  May  1983. 

R.  Perez  and  A.F.  Grandt,  "Initiation,  Growth,  and  Coalescence  of 
Small  Fatigue  Cracks  in  Polymethylmethacrylate,"  planned  for 
submission  to  5th  International  Congress  on  Experimental  Mechanics. 


R.  Perez,  "Cyclic  Growth  and  Coalescence  of  Multiple  Fatigue  Cracks 
Located  at  Notches,"  M.S.  Thesis,  Purdue  University,  expected 
May  1983. 
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Figure  1  Schematic  of  multiple  cracks  located  along  bore  of  a  hole 
In  a  large  plate.  Hole  diameter  Is  D,  plate  thickness  is  t 
and  remote  tensile  stress  is  applied  In  z  direction. 


Figure  2  Schematic  representation  of  semielliptlcal ,  corner,  and 
through  the  thickness  crack  located  along  the  bore  of  a 
hole  with  diameter  0  In  a  plate  with  thickness  t. 


Symmetric  Corner  Crack 


Configuration 


Configuration 
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Figure  3  Schematic  drawing  of  crack  plane  showing  location  of 
single  and  symmetric  corner  cracks  located  at  bore  of 
hole. 
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^  Hole  Bore 


INTERACTION  FACTORS 


Figure  4  Summary  of  the  effect  of  crack  spacing  tsep/a  and  crack 
shape  a/c  on  stress  Intensity  factor  Interaction  at  hole 
bore  crack  location  ($  *  90°). 
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VFree  surface 


INTERACTION  FACTORS 


Figure  5  Summary  of  the  effect  of  crack  spacing  tsep/a  and  crack 
shape  a/c  on  stress  Intensity  factor  interaction  at  front 
surface  location  ($  ■  0)  for  corner  cracked  hole. 
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Figure  7b  Comparison  of  Predicted  and  Actual  Crack  Growth 
During  Transition  Period  for  Plate  Specimen  T2 
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coalescing  corner 
in  cyclic  tension 


Figure  9  Comparison  of  predicted  and  experimental  fatigue  crack  growth 
for  two  coalescing  corner  cracks  located  at  a  hole  in  a  PMMA 
plate  loaded  in  remote  tension. 


c.  Crack  dimensions  after  coalescence 

Fig.  10  Specimen  geometry  showing  location  of  the  camera  and  schematic 
views  of  the  crack  plane. 
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Figure  11  continued 


CRACK  LENGTH  (  INCHES  ) 


Fig.  12 


THOUSANDS  OF  CYCLES  (  N  ) 

Crack  growth  before  coalescence  as  a  function  of  elapsed  cycles 
for  nonsymmetrlc  out-of-plane  crack  test  T-2. 
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CRACK  LENGTH  (  INCHES  ) 


CRACK  LENGTH  ( INCHES  ) 


THOUSANDS  OF  CYCLES  (  N  ) 

Fig.  13b  Crack  growth  after  coalescence  as  a  function  of  elapsed 
cycles  for  symmetric  In-plane  crack  test  T-3. 
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Fig. 14a  Crack  growth  before  coalescence  as  a  function  of  elapsed  cycles 

for  nonsymmetric  out-of-plane  crack  test  T-4. 
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Fig.  14b  Crack  growth  after  encounter  and  coalescence  as  a  function  of 
elapsed  cycles  for  nonsymmetrlc  out-of-plane  crack  test  T-4. 
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THOUSANDS  OF  CYCLES  (  N  ) 

Fig. 15  Crack  growth  before  coalescence  as  a  function  of  elapsed 
cycles  for  symmetric  out  of  plane  crack  test  T-5. 


Fig.  16  Crack  configuration  after  the  out  of  plane  overlap 
occured  before  coalescence. 
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DOUBLE  NOTCH  LOAD  CONTROL  LCF  SPECIMEN 


Comparison  of  predicted  and  experimental  growth  and  coalescence  of 
two  surface  cracks  at  semicircular  notch  in  Waspaloy  specimen 
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Figure  20  Comparison  of  predicted  and  experimental  growth  and  coalescence  of 
two  surface  cracks  at  semicircular  notch  in  Waspaloy  specimen- 
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Figure  21  Comparison  of  predicted  and  experimental  growth  and  coalescence  of 
two  surface  cracks  at  semicircular  notch  in  Waspaloy  specimen 


SPECIMEN  17 

INTERACTING  SURFACE  CRACKS 


Figure  22  Comparison  of  predicted  and  experimental  growth  and  coalescence  of 
two  surface  cracks  at  semicircular  notch  in  Waspaloy  specimen 


larison  of  predicted  and  experimental  growth  and  coalescence  of 
surface  cracks  at  semicircular  notch  in  Waspaloy  specimen. 


SPECIMEN -CAMERA  ARRANGEMENT 
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b.  CRACK  INITIATION  AT  V-NOTCH 


c.  INITIATION  SITES  UNITED  INTO 
UNIFORM  CRACK  FRONT 


Figure  24  Schematic  view  of  fatigue  crack  Initiation  at 
V-notch  and  coalescence  into  a  single  flaw 
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Figure  25a  Growth  of  small  fatigue  cracks  at  V-notch 
in  PMMA  bend  specimen 
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Figure  25b  Growth  of  small  fatigue  cracks  at 
V-notch  in  PMMA  bend  specimen 
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Figure  25c  Growth  of  small  fatigue  cracks  at 
V-notch  in  PMMA  bend  specimen 
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Figure  25d  Growth  of  small  fatigue  cracks  at 
v-notch  in  PMMA  bend  specimen 
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Figure  26a  Growth  of  small  fatigue  cracks  at  V-notch 
in  PMMA  bend  specimen 
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Figure  26b  Growth  of  small  fatigue  cracks  at  V-notch  in 
PWA  bend  specimen 
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Figure  26c  Growth  of  small  fatigue  cracks  at  V-notch  in 
PMMA  bend  .-'cimen 
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Figure  26d  Growth  of  small  fatigue  cracks  at  V-notch  in 
PMMA  bend  specimen 
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Figure  28  Growth  of  small  through- thickness  crack  at 
V -notch  in  PMMA  bend  specimen 
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Stress  Intensity  Factors  for  Coalescing 
and  Single  Corner  Flaws  along  a  Hole  Bore  in  a  Plate 

by 
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and 
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ABSTRACT 

The  objective  of  this  paper  is  to  describe  the  effects  of  crack 
interaction  on  stress  intensity  factors  for  two  symmetric  coplanar  corner 
flaws  located  along  a  hole  bore.  This  numerical  analysis  employs  the 
Finite  Element-Alternating  Method  to  determine  Mode  I  stress  intensity 
factors  for  single  and  coalescing  corner  flaws.  Using  single  flaw  stress 
intensity  factors  as  a  reference,  analysis  of  crack  size  and  shape  effects 
on  Kj  for  coalescing  corner  flaws  indicates  the  stress  intensity  factor 
for  crack  points  along  the  hole  bore  increases  as  the  crack  tip  separation 
distance  decreases.  Interaction  effects  are  not  experienced  by  hole  bore 
crack  points  when  the  crack  tip  separation  distance  is  equal  to  or  greater 
than  half  of  the  largest  corner  flaw  dimension. 
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INTRODUCTION 

Naturally  occurring  fatigue  cracks  frequently  initiate  at  several 
independent  points  along  the  bore  of  fastener  holes.  These  individual 
flaws  then  grow  and  coalesce  into  a  single  dominant  crack  which  controls 
final  fracture.  Prediction  of  this  coalescence  phase  of  crack  growth 
requires  stress  intensity  factor  solutions  for  the  individual  cracks  prior 
to  their  "link  up"  into  a  single  flaw. 

It  is  known  that  the  stress  intensity  factor  for  a  crack  tip  is 
influenced  by  the  close  proximity  of  an  adjacent  crack  tip  [1-4].  Crack 
interaction  effects  on  Mode  I  stress  intensity  factors  (Kj)  have  been  studied 
by  several  authors  [1-4].  Various  two-dimensional  Kj  solutions  exist 
for  interaction  of  coplanar  through-the-thickness  cracks  [3,  4].  Kamei 
and  Yokoburi  [3],  for  example,  present  a  two-dimensional  stress  intensity 
factor  solution  for  two  asymmetric  through-the-thickness  cracks  in  an 
infinite  sheet,  while  Benthem  and  Koiter  [4]  give  a  two-dimensional  Kj 
solution  for  two  symmetric  edge  cracks.  Murakami  and  Nemat-Nassar  [1] 
present  an  approximate  three-dimensional  Kj  solution,  obtained  by  the  Body 
Force  Method,  for  coplanar  surface  flaws  in  an  unbounded  solid.  None  of 
these  solutions,  however,  consider  coalescence  of  coplanar  flaws  located 
along  the  bore  of  a  hole  in  a  finite  thickness  plate.  Since  multiple 
fastener  hole  cracks  often  occur  in  service  [5,  6],  the  present  paper 
examines  crack  coalescence  along  a  hole  bore. 

The  objective  of  this  paper  is  to  describe  stress  intensity  factor 
results  for  two  adjacent  cracks  located  along  the  bore  of  a  hole  in  a 
finite  thickness  plate.  Figure  1  shows  typical  plate  cross  sections  for 
the  symmetric  and  single  crack  configurations  considered.  Here  a  is 
the  crack  dimension  measured  along  the  hole  bore;  c  is  measured  along  the 
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free  surface;  T  is  the  plate  thickness  and  D  is  the  hole  diameter.  The 
plate  is  loaded  with  a  remote  stress  a  applied  perpendicular  to  the  crack 
plane.  The  Finite  Element  -  Alternating  Method  (FEAM)  [7-13]  is  used  to 
compute  Kj  for  symmetric  (coalescing)  and  single  corner  cracks  along  the 
hole  bore. 

Results  are  given  in  both  tabular  and  graphical  form  for  crack  sizes 
in  the  range  0.25  <_  a/T  <  .5  for  crack  shapes  a/c  =  1.11,  1.5,  2.0  and 
3.0.  (As  described  in  earlier  work  [7],  numerical  programming  difficulties 
prevent  analysis  of  quarter  circular  cracks,  a/c  =  1.0,  with  the  present 
computer  codes.)  Crack  shape  and  size  effects  on  Kj  variation  around  the 
crack  perimeter  are  presented  for  coalescing  corner  flaws.  Two  and  three- 
dimensional  stress  intensity  factor  solutions  fcr  other  coalescing  crack 
configurations  are  compared  with  the  stress  intensity  factors  generated 
for  the  corner  cracked  hole  geometry. 
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NUMERICAL  APPROACH 

Stress  intensity  factors  for  coalescing  cracks  along  a  hole  bore  are 
obtained  by  the  Finite  Element  -  Alternating  Method  (FEAM).  The  present 
parametric  study  makes  use  of  computer  codes  developed  by  Smith  and 
Kullgren  [7].  The  original  codes  were  modified  here  to  treat  the  coalescing 
crack  problem. 

The  FEAM  calculates  Kj  based  on  approximate  surface  crack  boundary 
conditions.  An  iterative  superposition  of  two  solutions,  one  a  3-D 
finite  element  solution  for  an  uncracked  body  under  prescribed  surface 
loadings,  and  the  second  a  stress  solution  for  a  flat  elliptical  crack  in 
an  infinite  body  with  nonuniform  prescribed  surface  pressure, produce  these 
approximate  surface  flaw  boundary  conditions.  Earlier  applications  of  the 
FEAM  to  various  other  cracked  hole  problems  are  described  in  Ref.  8-12. 

The  crack  coalescence  problem  studied  here  requires  modification  of 
the  finite  element  mesh  and  the  elliptical  crack  stress  solution  of  [7] 
to  produce  the  symmetric  crack  finite  element  model  used  in  conjunction 
with  the  coalescence  option  of  the  FEAM  computer  code. 

Single  Crack  Model 

The  single  crack  model  is  a  semi-circular  plate  with  a  hole.  As  in 
Fig.  2,  this  model  has  a  hole  radius  R.  In  order  to  minimize  boundary 
effects  in  the  region  of  the  crack,  the  same  outer  plate  radius  of  12  R 
(R  =  hole  radius)  used  by  Smith  and  Kullgren  [7]  is  used.  The  single 
crack  model  employs  112  20-Node  isoparametric  finite  elements.  The  hole 
diameter  to  thickness  ratio  (D/T)  is  0.5.  The  modulus  of  rigidity  for 
the  plate  is  G  s  12  x  10®  psi  with  a  Poisson's  ratio  of  v  *  0.3. 


Symmetric  interacting  cracks  were  studied  by  imposing  a  line  of 
symmetry  along  the  x  =  0.5  face  of  the  original  FEAM  model  with  D/T  =  1.0 
as  shown  in  Fig.  2.  The  hole  radius,  outer  radius,  and  element  number  are 
equivalent  to  the  single  crack  model.  In  order  to  account  for  the  line 
of  symmetry,  the  symmetric  crack  model  has  a  D/T  =1.0  rather  than  D/T 
=  0.5  for  the  single  crack.  As  shown  in  Fig.  3,  the  plane  of  symmetry 
imposed  along  the  front  surface  of  the  plate,  however,  effectively  doubles 
the  symmetric  crack  model  plate  thickness.  Thus,  the  symmetric  crack 
model  D/T  equals  that  of  the  single  crack  model  (D/T  *  0.5).  Figure  3 
shows  how  an  infinite  plate  with  two  symmetric  cracks  along  the  hole 
bore  is  simulated  when  a  corner  flaw  is  applied  to  the  FEAM  symmetric 
crack  model. 


DISCUSSION  OF  RESULTS 


Stress  intensity  factors  were  generated  for  crack  shapes  a/c  =  1.11, 

1.5,  2.0,  and  3.0  by  the  FEAM  analysis.  Each  crack  varied  in  size  from 
0.25  <_  a/T  <  0.5.  The  stress  intensity  factors  for  the  corner  flaw 
configurations  considered  are  listed  in  tables  1  and  2.  Table  1  lists 
the  stress  intensity  factor  variation  for  the  single  crack  in  15  degree 
increments  of  the  elliptic  angle  $  defined  in  Fig.  1.  Stress  intensity 
factor  results  for  the  symmetric  crack  configuration  are  listed  in  Table  2. 

These  results  provide  a  parametric  study  of  size  and  shape  effects 

on  the  stress  intensity  factor  variation  around  the  crack  perimeter. 

For  a  single  corner  crack  with  a/c  *  1.11,  Figure  4  shows  the  variation 

of  the  dimensionless  stress  intensity  factor  (Kj/o/DT)  versus  normalized 

elliptic  angle  (<f>/4>  as  the  crack  size  (a/T)  increases.  Note  that 
max 

<t>max  =  90  degrees.  Figure  5  shows  the  variation  of  Kj/o/IT  over  the  same 
a/T  range  for  two  symmetric  cracks  with  the  same  a/c  *  1.11. 

Recalling  that  4>  *  90°  corresponds  to  the  hole  bore  crack  location, 
notice  the  increased  magnitude  of  Kj/o/D  along  the  hole  bore  for  the 
symmetric  configuration  compared  to  the  same  location  for  the  single 
flaw  as  the  crack  size  is  increased.  Using  the  single  flaw  Kj/a/U  as 
a  reference,  the  influence  of  a/T  on  Kj/o/U  for  particular  flaw  points 
on  the  symmetric  model  crack  perimeter  is  shown  in  Figs.  6-8.  For  small 
crack  sizes  (0.25  <_  a/T  <_  0.3)  shown  in  Figs.  6-8,  the  single  and  coalescing 
configurations  differ  by  a  maximum  of  4  percent.  However,  Figures  6  and  7 
show  a  significant  increase  in  the  symmetric  crack  Kj/o/U  as  compared  to 
the  single  flaw  Kj/o/U  when  the  crack  size  approaches  coalescence  (a/T  *  0.5). 
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Following  Murakami  [1],  a  crack  interaction  factor  y  is  defined 
as  the  ratio  of  the  stress  intensity  factor  for  two  cracks  to  the  stress 
intensity  factor  for  a  single  crack.  Kamei  and  Yokoburi  [3]  have  given 
two-dimensional  stress  intensity  factor  solutions  for  two  asymmetrical 
through-the- thickness  cracks  in  an  infinite  elastic  sheet.  By  equating 
the  crack  lengths  as  shown  in  Figure  9a,  and  normalizing  the  crack  tip 
stress  intensity  factor  expression  by  a/irT,  the  interaction  factor  y 
for  two  symmetric  through-the-thickness  cracks  is  obtained.  The  Benthem 
and  Koiter  [4]  two-dimensional  stress  intensity  factor  solution  for  two 
symmetric  edge  cracks  (see  Figure  9b)  may  be  used  to  compute  the  edge 
crack  interaction  factor  shown  in  Fig.  10.  Here,  the  double  edge-crack 
solution  was  divided  by  the  stress  intensity  factor  given  by  Harris  [14] 
for  a  single  edge-cracked  sheet  which  doesn't  bend. 

Three-dimensional  crack  point  interaction  factors  for  symmetric 

semi-circular  surface  flaws  (Fig.  9c)  are  presented  by  Murakami  and 

Nemat-Nassar  in  Ref.  [1],  These  coplanar  flaws  exist  in  an  unbounded 

solid.  The  present  FEAM  results  were  used  to  compute  three-dimensional 

interaction  factors  for  symmetric  corner  cracks  along  a  hole  bore  (Fig. 

9d).  Figure  10  compares  these  two-and- three  dimensional  crack  tip 

interaction  factors  as  a  function  of  dimensionless  crack  separation 
t 

distance  (— ) .  The  through-the-thickness  interaction  factor  of  Kamei 

fl 

[3]  represents  the  limiting  value  of  two  interacting  surface  cracks, 
while  the  Benthem  and  Koiter  [4]  edge-crack  factor  serves  as  a  limiting 
value  for  the  present  study  of  interacting  corner  cracks.  Due  to  a 
numerical  instability  in  the  FEAM  code,  values  of  crack  Interaction 
factors  at  coalescence  (tsep/a  ■  0)  are  not  obtainable. 


Notice  the  close  proximity  of  the  present  FEAM  results  for  corner 
cracked  holes  to  the  Murakami  [1]  interaction  factor  for  surface  cracks. 
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As  the  crack  separation  t$ep  exceeds  approximately  half  of  the  crack 
length  dimension  a,  (t$ep/a  >.0.5),  the  FEAM  crack  interaction  factor 
stabilize  at  unity,  indicating  that  the  crack  tip  located  along  the  hole 
bore  is  no  longer  influenced  by  the  adjacent  crack.  For  a  given  tsep/a, 
the  maximum  difference  in  all  hole  bore  interaction  factors  presented  in 
Figure  10  is  10.7  percent.  Figure  10  indicates  that  crack  size  (a/T) 
dominates  the  increase  in  Kj  due  to  crack  interaction  at  the  hole  bore, 
while  crack  shape  (a/c)  has  only  a  slight  influence.  The  effect  of 
crack  shape  on  the  hole  bore  intei action  factors  is  shown  on  an  expanded 
scale  in  Figure  11. 

Figure  12  presents  the  interaction  factor  at  the  plate  surface 
(♦  =  0).  The  factors  decrease  slightly  below  unity  in  the  range  of 
0  <  <  0.5.  Figure  12  further  indicates  that  the  crack  is  free  of 

influence  from  the  adjacent  crack  when  the  hole  bore  separation  distance 
is  greater  than  half  the  crack  length  a.  As  the  separation  distance 
decreases,  the  decrease  in  Kj  at  the  plate  surface  is  considerably  less 
than  the  increase  in  Kj  at  the  hole  bore  location  (see  Figure  10).  Thus, 
the  effect  of  crack  coalescence  is  fairly  localized,  serving  to  increase 
Kj  for  crack  points  in  the  vicinity  of  the  hole  bore.  The  free  surface 
Kj,  however,  is  only  slightly  affected  as  crack  separation  approaches 
zero.  Table  3  lists  the  hole  bore  and  free  surface  interaction  factors 
for  all  aspect  ratios  considered.  The  decrease  in  K  for  the  double 
cracked  case  may  be  caused  by  elimination  of  through- the -thickness 
bending  for  the  symmetric  crack  configuration. 
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SUMMARY  AND  CONCLUSIONS 


Several  conclusions  can  be  made  about  stress  intensity  factors  for 
two  symmetric  coalescing  corner  cracks.  First,  as  two  symmetric  corner 
cracks  approach  coalescence  (a/T  =  0.5),  the  greatest  increase  in  stress 
intensity  factor  occurs  at  the  crack  point  along  the  hole  bore.  The  crack 
point  at  the  free  surface  does  not  experience  a  significant  stress  intensity 
factor  increase,  but  shows  a  slight  decrease  in  Kj  as  the  crack  tips 
approach  along  the  hole  bore. 

For  crack  points  along  the  hole  bore,  crack  interaction  (stress 
intensity  factor  increase)  does  not  occur  until  the  crack  separation  distance 
is  equal  to  or  less  than  half  of  the  crack  length  a.  The  localized  effect 
of  coalescence  is  found  to  depend  strongly  on  crack  size  (a/T)  and  only 
weakly  on  crack  shape  (a/c).  The  crack  shape  dependence  indicates  that  the 
deeper  cracks  (smaller  a/c  ratios)  have  a  slightly  larger  K  (maximum  10.7% 
difference)  at  the  hole  bore  location  than  more  shallow  flaws  (larger  a/c) 
with  the  same  crack  length  a. 
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Table  1  Dimensionless  stress  intensity  factors  for  single  corner  crack 
at  open  hole  in  a  large  plate  loaded  in  remote  tension.  Results 
given  as  a  function  of  crack  shape  a/c, crack  size  a/T,  and  posi¬ 
tion  along  crack  perimeter  defined  by  parametric  angle  t 
($  *  0°  at  front  surface  and  =  90°  at  hole  bore). 


a/c  a/T 


K  j/a^“ 


.3000 

.3500 

.4000 

.4500 

.4750 

.4875 

.4938 


.2500 

.3000 

.3500 

.4000 

.4500 

.4750 

.4875 

.4938 


.2500 

.3000 

.3500 

.4000 

.4500 

.4750 

.4875 

.4938 


30°  45°  60°  75°  ♦  *  90" 


I 


.156 

.238 

.299 

1.388 


1.479 

1.496 

1.519 

1.540 


1.184 
1.270 
1.317 
1 .394 
1.448 
1.482 
1.519 
1.531 


.325 

.404 

1.426 

.440 

.456 


.137 
.183 
.234 
.285 
1.347 
1.371  1.422 

1.378  1.428 


flBSIE 


.398  1.571 

1.422  1.592 


1.542  1.773 

1.571  1.794 


1.388  1.433 


.150 

.268 

.308 

1.400 

1.435 

1.458 

1.480 

1.486 


1.313 

1.372 

1.410 

1.432 

1.444 

1.447 


.247 

.284 

1.315 

1.351 

1.386 

1.407 

1.414 

1.413 


.337 

1.385 

1.411 

1.427 

1.437 

1.439 


.302 

.330 

.365 

1.400 

1.432 

1.451 

1.460 

1.460 


1.813 


1.607 

1.610 


.460 

1.523 

.544 


1.577 


.532 

1.542 

.562 

.563 


.457 

.535 

.591 

1.676 


1.695 

1.712 

1.741 

1.743 


.596 

1.615 

1.646 

1.651 


1.665 

1.765 

1.831 

1.933 

1.952 

1.968 


1.840  2.007 

1.844  2.014 


1.676 

1.780 

1.799 

1.819 

1.859 

1.863 


.645 

.665 

1.705 

1.711 


Table  1  Continued 
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Table  2  Dimensionless  stress  intensity  factors  for  symmetric  corner 
cracks  located  at  intersection  of  hole  bore  with  front  and 
back  surfaces  of  a  plate  loaded  in  remote  tension.  Results 
given  as  crack  tips  approach  coalescence  (a/T  *  0.5)  in  term's 
of  crack  shape  a/c,  crack  size  a/T,  and  position  along  crack 
perimeter  defined  by  parametric  angle  $  =  90°  at  hole  bore) 


a/c  a/T 


Symmetric  Cracks  Kj/o  /S~ 


1.11  .2500 

.3000 
.3500 
.4000 
.4500 
.4750 
.4875 
.4938 


.2500 

.3000 

.3500 

.4000 

.4500 

.4750 

.4875 

.4938 


.121 

.185 

.210 

.276 

1.349 

.350 

.364 

.377 


.171 

.237 

.265 

1.321 

1.357 

1.384 

1.417 

1.426 


.145 

.201 

.224 

.283 

1.347 

1.356 

.356 

.347 


.219 

.273 

.301 

.352 

.385 

1.409 

.425 

.406 


30°  45* 


.182 

.229 

.262 

1.317 

1.381 

1.415 

1.425 

1.424 


1.644 


.254  1.317 

.296  1.358 

.332  1.410 

.382  1.470 

1.426 
1.458 
1.474 
1.464 


<T> 

O 

e 

75° 

1.422 

1.555 

1.479 

1.627 

1.560 

1.733 

1.643 

1.841 

1.765 

2.001 

1.872 

2.137 

1.981 

2.300 

2.127 

2.534 

1.457 

1.450 

1.526 

1.523 

1.617 

1.601 

1.715 

1.709 

1.854 

1.784 

1.948 

1.872 

2.081 

2.010 

2.302 

1.595 

1.680 

1.799 

1.925 

2.107 

2.250 

2.441 

2.722 


1.545 

1.645 

1.756 

1.912 

2.014 

2.173 

2.444 
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Table  3  Dimensionless  Interaction  factor,  y,  defined  as  ratio  of  stress 
Intensity  factor  for  symmetric  double  crack  geometry  divided  by 
single  comer  crack  result.  Interaction  factors  given  at  hole 
bore  ($  *  90°)  and  free  surface  (♦  3  0°)  locations  as  functions 
of  dimensionless  separation  along  hole  bore  U$ep/a)  and  crack 
shape  a/c. . 


INTERACTION  FACTORS  y 

HOLE  BORE  (* 

=  90°) 

FREE  SURFACE  (4*0°) 

a/c 

1.11 

1.5 

2.0 

3.0 

1.11 

1.5 

2.0 

3.0 

t  /a 
sep' 

.025 

1.352 

1.312 

1.278 

1.228 

..894 

.931 

.957 

.988 

.051 

1.216 

1.169 

1.136 

1.092 

.898 

.932 

.956 

.984 

.105 

1.143 

1.107 

1.084 

1.053 

.902 

.934 

.958 

.987 

.222 

.  1.079 

1.062 

1.051 

1.035 

.912 

.937 

.965 

.987 

.500 

.996 

.987 

.986 

.982 

.919 

.948 

.964 

.990 

.857 

.983 

.982 

.987 

.995 

.932 

.961 

.979 

1.009 

1.333 

.952 

.955 

.958 

.963 

.957 

.974 

.999 

1.020 

2.000 

.958 

.964 

.979 

.980 

.969 

.989 

1.033 

1.041 

Symmetric  Corner  Crack 
Configuration 


— 1c  h~ 


Single  Corner  Crack 
Configuration 

I- —  D  H 


Figure  1  Schematic  drawing  of  crack  plane  showing  location  of  single 
and  symmetric  corner  cracks  located  at  bore  of  hole. 


ill  il’A * 


a/T 

Figure  6  Comparison  of  dimensionless  stress  Intensity  factors  for 
single  and  symmetric  cracks  with  aspect  ratio  a/c  =  1.11 
showing  magnification  In  stress  Intensity  factor  at  hole 
bore  location  (<t>  *  90°)  as  crack  size  a/T  increases. 


K/aVD 


1 


3.5 


2.9 


2.3 


fl/C=l.ll  PHI=H5 


1.7 


S 


symmetric 

S 


_i _ i _ i _ i _ i 

.30  .35  .10  .15  .50 


a/T 

Figure  7  Comparison  of  dimensionless  stress  intensity  factors  for 
single  and  symmetric  cracks  with  aspect  ratio  a/c  »  1.11 
at  crack  perimeter  location  41  =  45°. 


K/on/O 


3-Si 


2.9 


2.3 


fl/C=l.ll 


i-dh 


singl 


.25  .30 
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Figure  8  Comparison  of  dimensionless  stress  intensity  factors  for 
single  and  symmetric  cracks  with  aspect  ratio  a/c  *  1.11 
showing  slight  decrease  in  symmetric  crack  stress  Intensity 
factor  at  surface  location  (<f>  *  0°)  as  crack  size  a/T  Increases. 


Figure  9  Interacting  crack  configurations  represented  by:  a.  Kamel 

(through- the- thickness  center  cracks  in  a  large  sheet),  b.  Koiter 
and  Benthem  (through-the-thickness  edge  cracks),  c.  Murakami  and 
Nemet-Nasser  (surface  cracks  In  a  semi-infinite  solid),  and  d.  present 
study  (symmetric  comer  cracks  located  along  hole  bore). 
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Figure  10  Comparison  of  effect  of  crack  spacing  t$ep/a  on  stress  Inten¬ 
sity  factor  for  various  two-  and  three-dimensional  crack  con¬ 
figurations.  The  Interaction  factor  y  Is  defined  as  ratio  of 
double  to  single  crack  stress  intensity  factors. 
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Figure  11  Summary  of  the  effect  of  crack  spacing  t$ep/a  and  crack  shape 
a/c  on  stress  intensity  factor  interaction  at  hole  bore  crack 
location  (*  =  90°). 
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Figure  12  Summary  of  the  effect  of  crack  spacing  t$ep/a  and  crack  shape 
a/c  on  stress  intensity  factor  interaction  at  front  surface 
location  ($  *  0)  for  comer  cracked  hole. 


